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Abstract 

The problem of enumerating spanning trees on graphs and lattices is considered. We 
obtain bounds on the number of spanning trees N$t and establish inequalities relating the 
numbers of spanning trees of different graphs or lattices. A general formulation is presented 
for the enumeration of spanning trees on lattices in d > 2 dimensions, and is applied to the 
hypercubic, body-centered cubic, face-centered cubic, and specific planar lattices including 
the kagome, diced, 4-8-8 (bathroom-tile), Union Jack, and 3-12-12 lattices. This leads to 
closed-form expressions for N$t for these lattices of finite sizes. We prove a theorem con- 
cerning the classes of graphs and lattices C with the property that N S t ~ exp(n2£) as the 
number of vertices n —>■ oo, where zc is a finite nonzero constant. This includes the bulk 
limit of lattices in any spatial dimension, and also sections of lattices whose lengths in some 
dimensions go to infinity while others are finite. We evaluate zc exactly for the lattices we 
considered, and discuss the dependence of zc on d and the lattice coordination number. We 
also establish a relation connecting zc to the free energy of the critical Ising model for planar 
lattices C 



1 Introduction 



The enumeration of spanning trees on a graph or lattice is a problem of long-standing 
interest in mathematics |J - [|J and physics 0-0- Let G = (V,E) denote a connected 
graph (without loops) with vertex (site) and edge (bond) sets V and E. Let n = v(G) — \V\ 
be the number of vertices and e(G) = \E\ the number of edges. A spanning subgraph G' is 
a subgraph of G with v(G') = \V\, and a tree is a connected subgraph with no circuits. A 
spanning tree is a spanning subgraph of G that is a tree (thus e(G') = n — 1). The degree 
of a vertex is the number of edges attached to it (often denoted coordination number or 
valence). A k- regular graph is a graph with the property that each of its vertices has the 
same degree k. For these and further related definitions see, e.g., [Q - 0. 

Denote the number of distinct spanning trees of a graph G by Nst{G). This number 
can be enumerated in terms of standard graph-theoretic quantities. Two methods for doing 
this will be used here: (i) via the Laplacian matrix [I], |2|, |], and (ii) as a special value of 
the Tutte polynomial || - [[10]. For the first, we recall the definition that two vertices are 
adjacent if they are connected by an edge. The adjacency matrix A of G is an n x n matrix 
whose elements are 

Aij = 1, if sites i and j are adjacent 

= 0, otherwise, (1.1) 

and the degree matrix A of G is an n x n diagonal matrix with elements 

= K{5ij , (1-2) 

where is the degree of site i, and Sij the Kronecker delta function. Define the Laplacian 
matrix 

Q = A - A . (1.3) 

Here and throughout this paper, we use boldface to denote matrices of size n x n. 

Since the sum of the elements in each row (or column) of Q vanishes, one of the eigenvalues 
of Q is zero. Denote the remaining n — 1 eigenvalues by Ai, A n _i. Then two basic theorems 
in graph theory state that 0, f| 



N ST {G) = Any cofactor of Q (1.4) 

(1.5) 



1 n—1 
n J- J- 



n . 
i=i 



An elementary proof of the equivalence of (|1.4j ) and (|1.5| ) can be found in |pT[ . The Lapla- 



cian matrix Q is also known in the literature as the Kirchhoff matrix, or simply the tree 
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matrix, which arose in the analysis of electric circuits || . The enumeration of N$t has very 
recently been considered in JTl|] for finite hypercubic lattices and the square net embedded 



on nonorientable surfaces. 

A second way to calculate Nst(G) is as the special value 



Nst(G) — T(G, 1, 1) (1.6) 



of the Tutte polynomial of the graph G p|j - fTO 



T(G,x,y)=Y,(x-V k(G ' ) - kiG) (y-V c(G ' ) i (1-7) 

G'CG 

where k(G) is the number of connected components and c(G) the number of independent 
circuits in G, and the summation is over all spanning subgraphs G' of G. Here, we have 
k(G) = 1 for connected graphs G, and it is clear that ( |1.7|) leads to ( p..6|) , since in the limit of 
x, y — > 1 the only contributing terms in (|1.7|) are those of k{G') = k{G) = 1 and c(G') = 0, 
namely, the spanning trees. 

In physics one often deals with lattices. A lattice is regular if all sites are equivalent [[12 1. 



For a wide class of graphs, including lattices which may or may not be regular and strips of 
lattices, the number of spanning trees N$t has the asymptotic exponential growth 

Nst(G) ~ exp(nz{ G }) as n — > oo . (1.8) 

Thus Z{g} provides a natural measure of the rate of growth, and is evaluated via 

Zs G } = lim -\nN ST (G) . (1.9) 

n^oo n 

where {G} denotes the formal n — * oo limit of a graph of type G. For lattices £ this is 
known as the bulk, or the thermodynamic, limit, and we denote z\g} by zc- Closed- form 
expressions for zc have been obtained for the square, honeycomb, and triangular lattices 
0, [7] . Exact results have also been obtained for strips of regular lattices of finite widths and 
infinite length [O - [ITH]. 



In the present work we consider spanning trees on general graphs G as well as in the 
limit of x, y — >• 1 lattices £ in d > 2 dimensions. Specifically, we shall 

(i) derive an exact relation between z^ and zc* for planar lattices £ and £* which are 
mutually dual. 

(ii) obtain bounds on Nst(G) and z^q\, 

(iii) establish the exponential growth Nst(G) ~ exp(nz{c}) for a wide class of graphs, 
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(iv) present a general formulation for the enumeration of Nst{£) and zc, 

(v) enumerate N$t{£) and evaluate z^ exactly for a number of lattices in d > 2, 

(vi) analyze the dependence of zc on the spatial dimensionality and the coordination 
number of the lattice C, including deriving an asymptotic expansion for d- dimensional hy- 
percubic lattices, and 

(vii) establish a relation connecting zc to the free energy of the critical Ising model for 
planar lattices L. 



Before proceeding, it is useful to review here the close connection of spanning trees with 
the Potts model in statistical mechanics The partition function of a 9- state Potts model at 



temperature T = 1//3 on G is |T§, [17 



Z(G, q ,v)=Y, e ~ PH 
R} 



1.10) 



where —f3Ti = K Yluj) $0*0$ with (ij) ranging over pairs of adjacent vertices in G. The 
summation is taken over Oi = 1, 2, q and % = 1, 2, n. The partition function fll.lOp can 
be written as (ITS 



Z(G,q,v)= J2l HG ' 



G'CG 



where the summation is over all spanning subgraphs of G and 



1 . 



1.11) 



1.12) 



The expression Ql.ll ) shows that Z(G,q,v) is a polynomial in q and v and enables one to 
generalize q from positive integers to real and complex values. The Potts model partition 
function on a graph G is related to the Tutte polynomial according to 



with 



so 



Z(G, q, v) = {x- l) k ^ G \y - 1)"T(G, x, y) 



x = 1 + q/v , y — v + 1 



1.13) 



1.14) 



q=(l-x)(l-y) . (1.15) 

Thus, Nst{G) can also be calculated in terms of the Potts partition function Z(G, q, v). The 
values x = 1 and y = 1 in the evaluation ( |1.6| ) correspond to the limits q — > 0, and v — > 
with lim { g^}_ (g/t;) = 0. 
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For a planar graph G and its dual graph G*, the Tutte polynomial possesses the duality 
relation g |, 

T(G,x,y)=T(G*,y,x). (1.16) 

An immediate consequence of ( |1.6| ) and ( |1.16| ) is that the number of spanning trees is the 
same for a planar graph G and its dual G* (see, for example, [Q]), namely, 

Nst{G) = Nst{G*) . (1.17) 

Let n* be the number of sites of G*, given by the Euler relation 

n* = \E\-n+l. (1.18) 

For planar lattices £ and its dual £*, it is convenient to introduce the ratio in the bulk limit 

u c = lim (n*/n) (1.19) 

n— >oo 

satisfying 

= 1 . (1.20) 

Using Q1.19Q and (|1.20|) , we obtain the relation 

zc* = z c /vc (1.21) 

relating 2£ and zc*. As examples, it is readily verified that we have 

v hc = 1/2, v kag = 1, z/ 4 _ 8 „ 8 = 1/2 , (1.22) 

where the subscripts denote the honeycomb, kagome, and 4-8-8 lattices, respectively (for a 
detailed discussion on classifications of planar lattices see, for example, [12|, fHJ). Applying 



(|1.21|) to the respective duals, namely, the triangular, diced, and Union Jack lattices, we 



obtain the relations 

Ztri — 2 Zhc , Zdiced = Zkag , Zjjj = 2 Z^8-8- (1-23) 

For a regular planar lattice involving a tiling of the plane with only one type of polygon, 
the dual lattice is also regular. This includes the square which is self-dual, and honeycomb 
and triangular lattices which are mutually dual. In contrast, for a regular planar lattice 
involving a tiling of the plane with more than one type of polygon, the dual lattice is not 
regular since it does not have a uniform coordination number. For example, the kagome and 
4-8-8 lattices, which are regular, involve tilings with more than one type of polygon, and 
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their duals, the diced and Union Jack lattices, are not regular. For nonregular lattices it is 
convenient to introduce an effective coordination number n e ff defined as the average number 
of edges per vertex, 

K eff = lim (2\E\/n) . (1.24) 
Combining (|1.18|) and ( |1.24| ), we derive the relation 

K eff = 2(l + v c ) . (1.25) 



As examples, using (|1.20| ) this yields « e // = 4 and 6, respectively, for the diced and the 
Union Jack lattices. For K-regular graphs we have K e ff = k. 

2 General Bounds 

The determination of upper bounds on Nst{G) is a problem of considerable interest in graph 



theory. A general upper bound is [21 



n \n — 1 

For a /t-regular graph, this implies the upper bound 



i / 21 E\ \ n ~ 1 

N ST (G)<-(-^-) , (2.1) 



n\n — 1 

and hence 



NsAG)<-(^-) n ' 1 (2.2) 



z c <lnK. (2.3) 

More generally, we shall also deal with lattices that are not regular but for which one can 
define an effective coordination number, K e ff as in ( |1.24| ). For these, from ( |2.1| ), one has 



z c < In K,.ff. (2.4) 

A stronger upper bound for K-regular graphs with k > 3 is due to Mckay, Chung, and 
Yau who established rigorously that 

where 

K - [k{k - 2)]«/2-i " 1 ' 
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This leads to the upper bound 



z c < InC, 




K 



12/t 3 4k 4 



7 3 



+ 



30k 5 



31 



2k 6 \k 7 



) 



(2.7) 
(2.8) 



where we have carried out a large-K expansion. 

On the other hand, lower bounds on Nst{G) are more difficult to obtain. We have 
established an inequality which we state as a theorem: 

Theorem 2.1. 

Let G be a connected graph and let i and j be two nonadjacent vertices of G with degrees 
Ki and Kj, respectively. Let H be a graph obtained from G by adding an edge 
connecting z and j. Then we have 



where K = minj^j, Kj}. 

Proof: Since sites i and j are connected in G, the adding of the edge to any spanning 
tree T(G) on G forms a closed circuit on H. The closed circuit contains in addition to 
another edge £ incident at site j. The deletion of I then breaks the circuit, resulting in a 
spanning tree configuration T(H) on H. However, the spanning tree T(H) so constructed 
is not necessarily unique; the same T(H) may result from m different T(G). Since each 
T(G) is also a spanning tree of H, it follows that for the m spanning trees T(G) there exist 
m + 1 distinct spanning trees T(H). A moment's reflection shows that we have always 
m < k, with Tn — k arising when there is a single edge incident at site j in T(G) and a 
single edge at j in T(H). Since we have (kj + l)/ft 3 - > («j + l)/«t for ^« > fy, the 
proposition follows as a consequence. □ 

Corollary 2.1. 

Let G be a K-regular graph, and let G' be a graph derived from G by adding M edges one at 
a time such that each added edge terminates in at least one vertex whose degree is k. Then 



Remark: Corollary 2.1 is proved by applying Theorem 2.1 M times. Furthermore, for 
lattices z c and z c >, and M = an where a is a constant, ( |2.10|) implies the bound 




(2.9) 




(2.10) 




(2.11) 
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For example, by adding edges one at a time one can convert the honeycomb lattice first to 
the square and then to the triangular lattice. Corollary 2.1 then implies the inequalities 



4 

z sq > z hc + - In I - 



and 



(2.12) 



(2.13) 



z tri > z sq + hi 

Combining these bounds with the relation Ztri = 2zh c from ( |1.23| ), we obtain the lower bounds 

0.510 825 6... (2.14) 
: 0.366 984 5... (2.15) 



z sq > In 



2^ 



z hc > In ( 
z^i > 2 In (^) = 0.733 969 1... 



(2.16) 



Note that these bounds have been deduced without actually carrying out explicit calcula- 
tions. For comparison, the exact values are || |7| (see also Sec. 4 below) 



"sq 



71 



3 2 + 52 7 2 + 92 



1.166 243 6... 



(2.17) 



and 



Zhc 
Ztri 



3^ 

~2t7V 
2z hc = 1.615 329 7. 



1 1 1 

72 "TP + 132 



0.807 664 9. 



(2.18) 
(2.19) 



To measure the effectiveness of these lower bounds, let R{g} denote the ratio of the exact 
z {G} to the bounds ( gig) - (gl6|) . We have 



Rsq ^ 2.28, 



R 



he 



Rtri ~ 2.20 



(2.20) 



indicating that the bounds are very generous. 

It is also of interest to work out the implications of the lower bound ( 2.11 ) for strips of 
regular lattices. For example, for 2 x 00 square [O, KM and triangular [|L5] strips one has 



the results [O, 111, IT 



4(2xoo) = ^m(3 + V2 ) = 0.658 478 9... 



free 
^tri(2xoo) 



iln(I±M) 



0.962 423 6... 



(2.21) 
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where the superscript denotes free boundary conditions in the transverse direction. This 
is to be compared with ( |6.6.1|) and ( |6.6.2|) below for periodic boundary conditions in the 



transverse direction for which we have z 



sq(2xoo) 



0.881 373 5... and z t 



ri(2xoo) 



1.386 294 3... 



It is easily checked that these numbers obey ( 2.13Q . Furthermore, the inequality ( 2.13 ) now 
implies 



1 (\ 

Z tri(2xoo) > Z sq(2xoo) 2 



free 



Jree 



Ztri(2xoo) > ^sg(2xoo) + 2^ n ( 



0.802 319 9... 
1.025 214 5... 



so that the ratios of the exact values to the lower bounds (|2.22|) are 



R 



free 

tri(2xoo) 



1.20, 



Rfr 



tri(2xoo) 



1.35 . 



(2.22) 



(2.23) 



Another example is the <i-dimensional hypercubic lattice d with k = 2d for which 
Corollary 2.1 implies^] 



z Cd+1 >z Cd +\n l + (2d)- 1 



(2.24) 



Applying this lower bound to z sq using zu ne = 0, we obtain z sq > ln(3/2), which is not as 
strong as the bound ( |2.14| ). For the d = 3 simple cubic lattice, we have 



z sc > z sq + In 



1.389 387 1. 



so that 



R sc ~ 1.20 . 



(2.25) 
(2.26) 



Similar lower bounds can be deduced for hypercubic lattices with d > 4. 



3 Families of Graphs with Exponential Growth for N$t{G) 

In this section we prove a result concerning the class of families of graphs for which Nst{G) 
has the exponential asymptotic behavior ( |1.8| ). 

A family of graphs is recursive if it can be built up by sequential additions of a given 
subgraph. As an example, consider a strip of width jV 2 and length iV\ = m of some lattice 
such as a square lattice; this can be built up by starting with a column of N 2 — 1 squares and 
sequentially adding columns to elongate the strip in the N\ direction. A higher-dimensional 



1 Here we add a technical remark. As the inequality (2.24) is deduced by constructing Cd+i via adding 
edges to connect n x / d copies of Cd, a technical problem arises if the copies are disjoint to begin with for 
which Nst = by definition. The difficulty is resolved if one starts instead from copies of Cd which are 
connected by a single edge. 
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example is a rectangular tube of a lattice such as a simple cubic lattice with transverse size 
N 2 x N 3 and length N\ = m; this can be built up by starting with a single N 2 x N 3 section 
and sequentially adding m transverse sections and connecting them in an obvious manner 



to elongate the tube. We will need the following result from [14 



Lemma 3.1. 

Let G m be a recursive graph of length m subunits as described in the above. Then the 
Tutte polynomial has the form 

N a 

T(G m ,x,y) = ^2c GJ (x,y)(a Gij (x,y)) m , (3.1) 

where the explicit forms of a G j(x,y) and c G j(x,y) depend on the type of graphs G m . 

The proof of Lemma 3.1 can be found in Ul4| . (See (2.18) and (8.15) of ||14fl , where a G ,j was 
denoted by X G j.) Note that the class of recursive graphs is more general than the class of 
regular lattices. 



As discussed in [H , as m — > oo for a given (x, y), the term a G j with the maximal magni- 



tude will dominate the right-hand side of ( ^.1| ), provided that the corresponding coefficient 
c G j(x,y) does not vanish. We denote this term a G ^ )7nax . Using the relation between the 
Potts partition function and the Tutte polynomial ( |1.13| ), together with the definition of the 
(reduced) Potts model free energy 

f({G}, g, v) = lim n" 1 In Z(G, q, v) , (3.2) 

n^oo 

one observes that a nonanalyticity in / occurs when, as one changes (x, y) or equivalently 
{q,v), there is a crossover of the dominant term a G j. These changes determine the regions 
of analyticity (phases) of the free energy. We next proceed to our theorem. 



Theorem 3.1. 

Let G m denote a recursive graph of a lattice of length m in one spatial dimension with 
fixed (d — l)-dimensional "transverse" section of the size ^ x Af 3 x ■ • Nd- For 
N 2 , ■ ■ •, Na > 2, the number of spanning trees NsT{G m ) grows exponentially with n as in 
( |1.8[ ), thereby defining a nonzero finite constant Z{ G y. 



Proof. We shall carry out the proof for the case d = 2; the generalization to d > 3 is 
straightforward. The strategy of the proof is to use the structural result ( |3.1|) in Lemma 
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3.1 above. It is convenient to use some results from the Potts- Tutte correspondence (|1.13| ). 
From (|1.12| ) and (|1.14|) it follows that y = 1 corresponds, in terms of the correspondence 



( L13 ), to the infinite-temperature point K = v = 0. From the basic property that a spin 
model such as the Potts model is analytic at K = 0, i.e., has a Taylor series expansion in K 
(or v ) with a finite radius of convergence, it follows that, in the neighborhood of y = 1 for a 
given x, a single term 

0>G,max = OG.fli (3.3) 

will dominate, where R\ denotes the region in the (x, y) space corresponding to the 
paramagnetic phase for a given q in the (q, v) space. It was shown in |TJ| that the 
corresponding coefficient in T{G,x,y) is nonzero. Further, one knows that aG jmax (l, 1) > 1 
unless G is the tree or circuit graph or obvious modifications thereof, a fact which can be 
easily proved by assuming the contrary and deducing that N$t violates a lower bound on 
the number of spanning trees for strip graphs with width > 2. From (|1.6|) and ( |1.9|) , it 
follows that NsriGm) has the exponential asymptotic growth as in ( |1.8|) , so that there is a 
finite nonzero constant ^{g}- The method of the proof evidently provides a constructive 
way to calculate this quantity in terms of (|3.1|), or 



z {G} = t 1 ln[a Gd>max (l, 1)] , (3.4) 

where t is the number of vertices in a transverse section. It is straightforward to generalize 
these results to the case of a <i-dimensional tube graph with fixed (d — l)-dimensional 
transverse section. This completes the proof. □ 



To place this result in perspective, we give some examples of families of graphs for which 
the asymptotic behavior ( |1.8|) does not hold. For the tree graph T n and circuit graph C n , the 
Tutte polynomials are T(T n , x, y) = x n ~ l and T(C n , x, y) = y+J2s=i %s so ^ na ^ ^ST{T n ) = 1, 
independent of n, and NsriCn) = n. In both cases z^t} = z {c} — 0. It is for this reason 
that we restricted the width of the strips to A^2, N3, ■ ■ ■ > 2 in Theorem 3.1. An example of 
a family for which Nst grows more rapidly than an exponential is the complete graph K n , 
a graph with the property that each vertex is adjacent to every other vertex. In this case, 
one has N ST (K n ) = ri n ~ 2 0. Note that K n is not a recursive graph. We next proceed to 
the case of the bulk limit of a lattice. 

Theorem 3.2. 

Let C denote the bulk limit of a regular rf-dimensional lattice of size Ni x ... x Nj, with 
Nj — > 00, j = 1, d such that the ratios lim^oo Ni/Nj remain nonzero and finite. Then 
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the number of spanning trees on this lattice grows exponentially with n as in (|1.8|) , thereby 
defining a nonzero finite constant Zjr. 

Proof. A sketch of the proof is given here. The idea is to use Theorem 3.1 and observe that 
the property of exponential growth of N$t as ^1 — > 00 f° r fixed N 2 is independent of 
N 2 > 2. Furthermore, as discussed in the proof of Theorem 3.1, given the relation ( |1.6| ), it 
follows that Nst and z are, in statistical mechanics terminology, determined by Z(G,q,v) 
and the free energy / at the disorder point v = K = 0. Hence, we can take the limit 
N 2 —>■ 00, and since the exponential growth of Nst holds uniformly in N 2 , it also holds in 
this limit. This establishes the result for d = 2. It is straightforward to extend the proof to 
d > 3 by using Theorem 3.1 starting with a tube of the ci-dimensional lattice with fixed 
(d — l)-dimensional transverse cross section and one longitudinal direction that goes to 
infinity, and using again the fact that the exponential growth of Nst holds uniformly as 
one increases the (d — 1) transverse dimensions of the tube. □ 



We remark that Theorem 3.2 can also be established directly using the explicit expression 
( |4.12j ) of Nst{£) obtained in the next section. Another remark is that it is important that 
z is a disorder quantity. Because of this, one gets the same asymptotic behavior for z on 
a sequence of infinite-length tube graphs of progressively larger and larger transverse cross 
sections as one does by taking the limit of all Nj — > 00 with the ratios N/Nj fixed. In 
contrast, this would not be the case in dealing with a quantity connected with a divergent 
correlation length. For example, the free energy of the ferromagnetic Potts model on infinite- 
length, finite- width strips has a zero-temperature critical point at v — > 00 for any finite width, 
but has a quite different analytic structure if one takes both N\ and N 2 to infinity with N 2 /Ni 
fixed, namely, a non-analyticity (phase transition) at a finite temperature v > 1. 



4 Formulation for General Lattices 

The formulation of enumerating spanning trees for general lattices is given in this section. 
Consider a lattice £ of n sites in d spatial dimensions. We shall use ( |1.5[ ) to evaluate Nst 
and for simplicity assume periodic boundary conditions. Formulations for other boundary 
conditions can be similarly worked out (see, for example, [ TT | ) . 



To write down the Laplacian matrix Q in a form suitable for computing its eigenvalues, we 
make use of the fact that any lattice in d dimensions is decomposable into a hypercubic array 
of N% x N 2 x • • • x Nj unit cells, each containing v sites so that we have n = vNiN 2 ■ ■ ■ A^.Q 
2 If Nj = 2 for some j, then the two sites in the j-th direction are connected by double edges. But this 
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Specify the cells by the coordinate n = {ni, n 2 , • • • , n^}, where 7ij = 0, 1, 2, • • • , Ni — 1, and 
number the sites in a cell 1, 2, • • • , v. Let a(n, n') be the v x v cell vertex adjacency matrix 
describing the connectivity between the vertices of the unit cells n and n'. Namely, 

aij(n, n') = 1, if site i in cell n and site j in cell n' are adjacent 

= 0, otherwise. (4-1) 

Under the assumption of periodic boundary conditions, we have the translational symmetry 

a(n, n') = a(n — n'), (4-2) 

and we can therefore write a(n) = a(n 1; n 2 , • • • , n d ). 

The general formulation is best illustrated by considering an example. Here we consider 
the 4-8-8 (bathroom-tile) lattice shown in Fig. 0(a). This is a regular lattice which has the 
coordination number 3, and the unit cells are the obliquely oriented squares, with v = 4. 
The explicit forms for the a(ni,n 2 ) matrices depend on one's convention for labeling the 
vertices within a unit cell; we choose the labeling shown in Fig. |l](a). Then one has 




(a) (b) 

Figure 1: (a) The 4-8-8 (bathroom tile) lattice, (b) The 3-12-12 lattice. Sites within a unit cell are labeled 
as shown. 



a(0,0) 



/ 1 1 \ 
10 10 
10 1 
\1 1 0/ 



does not affect any of the ensuing discussions. 



a(l,0) 



:-i,o) 



/ 1 \ 





\ / 



12 



a(0,l) 



a T (0, 



( \ 

1 



\0 0/ 



(4.3) 



where a T is the transpose of a. Then the Laplacian matrix Q assumes the form 

Q = 3J 4 -a(0,0) ®I Nl ®I N2 -a{l,0)®R Nl ®I N2 + a(-l,0)®R T Nl < 
-a(0, 1) ® I Nl <g> R N2 - a(0, -1) <g> ® i$ 2 , 

where In is the N x N identity matrix, and Rn the N x N matrix 



l N 2 



(4.4) 



/ 1 
1 



R 



N 



\ 





1 



• o\ 

• 

• 1 
■ °/ 



(4.5) 



To determine the eigenvalues of Q, we make use of the fact that Rn is diagonalized by 
the similarity transformation SnRnS^ 1 generated by the matrix Sn with elements 



(S N ) nm = (S^y mn = AT-VV 2 ™^, m, n = 0, 1, • • • , N - 1, 
where * denotes the complex conjugate, yielding the eigenvalues 

\ n = e ,2 ™/ N , n = 0,l,---,JV-l. 
It follows that the similarity transformation generated by 

I u ® S Nl (g) S N2 



(4.6) 



(4.7) 



diagonalizes Q in the and A*2 subspaces. Then, using the fact that a determinant is equal 
to the product of its eigenvalues, we obtain from (|1.5|) the expression 



A 



4iViA^ 2 



ATi-l N 2 -l 

n n« 

fci=0 k 2 =0 



M 



2irki 27r/c2 



(h,k 2 )^ (0,0), (4.9) 



where A = AiA 2 A 3 = 64, and Ai = A 2 = A 3 = 4 are the nonzero eigenvalues of M(0, 0). Here, 
M is a 4 x 4 matrix and the notation det|M| denotes the determinant of M. Explicitly, we 
have 



M(0i,0 2 ) 



3J 4 - a(0, 0) - a(l, 0)e ifll - a(-l, 0)e" tei - a(0, l)e iy2 - a(0, -l)e 



-i0 2 
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/ 


3 


-1 -e iSl 


-1 


\ 




-1 


3 -1 


_ e ie 2 






_ e -iBi 


-1 3 


-1 




V 


-1 




3 


/ 



and hence 



det|M(0i,0 2 )| =4 



7 — 3 (cos 6i + cos 62) — cos #1 cos 1 



(4.10) 



Note that a change in labeling conventions would either interchange or negate 9\ and/or 62- 
These have no effect on the final expression since det|M| is invariant under these changes. 
From ( 4.10 ) we have 

£(£4-8-8) 



-ln2 + ~ / — ^ I ^ln 



4 

4 47T 



0.786 684(1). 



2tt 



7 — 3 (cos 0i + cos 62) — cos 6\ cos 6*2 



d# In 



7-3cos# + 4sin(fl/2)V5 — cos 9 



(4.11) 



The last two lines are obtained after carrying out one integration followed by a numerical 
integration of the remaining integral. 



The consideration of general L now proceeds in a similar fashion. In place of 
obtains 



Nst(C) 



A 



vN x ■ -N d 

where k = (ki, k 2 , ■ ■ ■ , k d ), = (0, 0, • • ■ , 0) 



Ni-l N d -1 

n 

fe 1= k d =0 



2irk c 



(k^O). 



D{e u ---,e d ) = det \M{e x ,---,e d )\ , 

A = 1, v = \ 

= Ai • -\ v -i, v > 1 



one 
(4.12) 

(4.13) 
(4.14) 



and Aj's are the v — 1 nonzero eigenvalues of the matrix M(0, •■, 0). Here, M is a v x v matrix 
defined by 

M(6 1 , ■■■,6 d ) = A u -J2 a(n)e in ' e , (4.15) 

n 

where A„ is the degree matrix (|1.2|) for a unit cell, and — (9 



(9i, -,9 d ). Note that the 
determinant of M is always real since the matrix M is hermitian. This leads to the result 



v 



Zc = z (C) = - / -±... / -^\nD 



2tt 



* d6 H 



2tt 



(4.16) 
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Eqs. Q4.12D and (|4.16| ) are our main results for general regular lattices, and (|4.12|) is suitable 
for enumerating N$t for lattices of finite sizes. 

It is also of interest to consider the case in which the size of the lattice is finite in £ < d 
dimensions and goes to infinity in the remaining d — i dimensions. For example, if one sets 
I = 1, then as N± increases from 1 to infinity, the resultant sequence of values of z can be 
regarded as a sort of "interpolation" between the infinite (d — l)-dimensional and infinite 
ci-dimensional lattices. Without the loss of generality we let Li,---L,£ be finite. Then from 
Q4.12j ) we have 

1 



iVi-l 



• x Ng x oo x ■ 
n d0£+i 



xoo 



fci=0 



Nt-1 

E 

kt=0 



2vr 



71 d(h 



D 



27rfci 



x 



iVi 



27cki 



i+i-, 



e, 



(4.17) 



This shows clearly the result of Theorem 3.1, that N$t grows exponentially like ( |1.8| ). For 
d = 2,£ = 1, i.e., the case of iV x oo strips of infinite length and finite width Ni = N, the 
integration can be carried out. In particular, for strips of square and triangular lattices, we 
obtain the expressions 



z sq(N xoo) 



Ztri(Nxoc) 



1 



N 
1 

N 



COSUk 



1 



1/2 



fc=0 
N-1 



E'H 3 - 



fc=0 



COStUfc 



1 — cosu;fc)(7 — cos a;*/ 



1/2 



(4.18) 



(4.19) 



where u>k = 2irk/N. The summations Q4.18 ) and ( [4.19|) can be explicitly carried out for a 
given N. We give the results in Sec. 6.6 below. 



5 Lattices in d > 3 Dimensions 

The formulation of the preceding section is now specialized to specific lattices in d > 3 
dimensions. 



5.1 d-Dimensional Hypercubic Lattices 

Consider first the (i-dimensional hypercubic lattice Cd consisting of n = N1N2 ■ ■ ■ Nd sites. 
Here, £2 and £3 are the square and simple cubic lattices. We have v = 1, Ai = k — 2d, and 
the adjacency matrices 

a(n) = 1, n = (±l,0,---,0),(0,±l,---,0),---,(0,0,---,±l) 

= 0, otherwise. (5.1.1) 
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Using the general expressions Q4.12 ) and ( |4.16| ), one obtains immediately 



N ST (C a 



)n— 1 



n 



n 

fei=0 



From (|5.1.2j ) we derive 
z{£ d ) = ln(2rf) + 



N d -1 

n 

k d =0 



w d6i 
„~2^ 



d 



cos 



2fci7T 



+ 



+cos Ur) 



(k^O) (5.1.2) 



d0d 
2tt 



In 



1 - -(cos 0i + ■ 



+ cos 9 d ) 



(5.1.3) 



The expression (|5.1.2 ) is the same as that reported in [nj|. 

As in the case of the 4-8-8 lattice, one of the integrations can be carried out analytically 
and the remaining (d — l)-fold integral done numerically. For d = 2, the quantity -z(£ 2 ) can 
be exactly evaluated || leading to the value given by ( [2.17 ) which we include in (|5.1.5| ) 
below. We have carried out the numerical integrations for d — 3,4 and present the results 
also in ( |5.1.5|) . For higher values of d the evaluation of the integral becomes less accurate; 
however one can expand the logarithm in ( |5.1.3| ) and carry out the integrations term by 
term. This leads to the large-ci expression 

7 45 269 



z{C d ) = ln(2d) - 



1 3 
4d + 16cP 



805 /l 
• ) 1 2d 1 ' \d 7 J. 



(5.1.4) 



32d 3 128# 384d 5 

Note that the expansion ( J5.1.4 ) agrees with the corresponding large-K expansion ( |2.8| ) of the 
Mckay-Fan-Yau upper bound to the order of after taking k = 2d. However, 1/d expan- 
sions are generally expected to be asymptotic in view of the slow growth of the coefficients 
(see, for example, P^J|). If one truncates the 1/d series to a given order and lets d — > oo, 
the approximation to z(C d ) becomes progressively more accurate. But for a fixed d, one 
does not necessarily obtain a more accurate approximation by including more terms in the 
calculation. In practice, however, we found that, for d = 4,5, and 6, the series evaluated 
to 0(d~ 4 ) and 0(d~ 5 ) gives essentially the same values, which we listed in ( p. 1.5 ), with an 
accuracy of 10~ 4 . Combining our results, we have 



z(C 5 ) 
z{£ 6 ) 
z{C d ) 



1.166 243 6 - • ■ 
1.674 148 1(1) 
2.000 0(5) 
2.243 
2.437 

ln(2d), d 



(numerical evaluation) 
(numerical evaluation) 
(series expansion) 
(series expansion) 



oo. 



(5.1.5) 



Our numerical result suggests that z(C$) may be exactly equal to 2. It is readily verified 
that values of z(£ d ) in ( |5.1.5|) are consistent with the inequality (|2.24|) . 
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5.2 (/-Dimensional Body-Centered Cubic Lattice 



For the usual three-dimensional body-centered cubic (bcc) lattice, a unit cell contains v — 2 
vertices located at (0,0,0) and (~, 5, ~) numbered 1 and 2, respectively. Then one has 
A 2 = 8J 2 and the adjacency matrices 



a(0,0,0) 
a(l,0,0) 



a(-l,0,0) 



1 

1 



a(0,l,0) = a(0,0, 1) = a(l,l,0) 
a(l,0,l) = a(0,l,l) = a(l,l,l) = 
a(0, -1, 0) = a(0, 0, -1) = a(-l, -1, 0) 
a(-l,0,-l) = a(0,-l,-l) = a(-l,-l,-l) 



1 






1 



leading to the matrix 



M{e 1 ,e 2 



-VIV2V3 



-(viv 2 v 3 )* 



where Vj = 1 + e 10 ', j = 1, 2, 3. The evaluation of the determinant yields 



64 



1 - cos 2 (#i/2) cos 2 (fl 2 /2) cos 2 (fl 3 /2) 



This leads to 
N ST {C bcc ) 



iVi-l N 2 -l N 3 -l 

n n n 

fe 1= k 2 =0 k 3 =0 



64 - 64 cos 2 



ki7T 



COS 



k 2 ir 



cos 



%7T 



(^1,^,^)^(0,0,0) 



(5.2.1) 



(5.2.2) 



(5.2.3) 



(5.2.4) 



where n = 2N\N 2 Ns, and A = 16 is the nonzero eigenvalue of M(0,0). The expression 
( |5.2.4|) enumerates Nst for finite bcc lattices. Using (|1.9| ), we obtain 



bcc ) 



31n2 



31n2 + 



1 r de l r do 2 r do. 



2tt 



2tt 



2tt 



In 



1 — cos 



COS 



COS 



2tt 



2tt 



2tt 



In ( 1 — cos 61 cos 6 2 cos 6-- 



03 

(5.2.5) 



These results can be generalized to the d-dimensional body-centered cubic lattice, which we 
shall denote bcc(d). This lattice has coordination number k = 2 d . For finite lattices we 
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obtain 



NsT(£bcc(d)) 



A 



2Nt ■ -N ti 



Wi-l iV d -l 

n n 

ki=0 k d =0 



J] cos 2 

3=1 



kfK 



(MO) (5.2.6) 



where A = 2 d+1 . Taking the bulk limit, we obtain 



d i n2+ r^i... r^i n (i_( COS ^)... (cos ^ 

J-tt 2tt 2tt V 



d In 2 



2^ £\ 



If (2£)\ 



(5.2.7) 



where we have expanded the logarithm and carried out the integrations term by term. 

For d = 2, the body-centered cubic lattice bcc(2) is just the square lattice, and it is 
readily seen that the two integrals in ( |5.1.3| ) and ( [5.2.7] ) are equal for d = 2. Interestingly, 
this also establishes the equality of the two series in (2.17) and ( |5.2.7| ) at d = 2. We have 
further evaluated z(Cbcc(d)) for d > 3 using both expressions in ( |5.2.7| ), and found that the 
series converges slowly, with good agreement between the two reached only after evaluating 
the series to 100 — 200 terms. For d = 3 and 4 the results are 



z(C bcc ) = 1.990 2(1), 
z(£ M4) ) = 2.732 3(1), 



d = 3 
d = A. 



(5.2.8) 



5.3 Face-Centered Cubic Lattice 



For the face-centered cubic (fee) lattice, a unit cell contains v = 4 vertices located at (0, 0, 0), 
1 1 , J n 1 ' ; L 0) numbered 1,2,3,4, respectively. One has A 4 = (12) J 4 and the 



f0 i i) (± -) ( 



1 1 

2' 2' ' 



adjacency matrices, 
I ° 

a(0,0,0) = 

\ 

T I 



— 1 



-1 
-1 







a(0,l,0) = a 1 (0,-1,0) 



a(l,l,0) = a T (-l, -1,0) 



/0 
1 


V 1 



-i \ 
-i 
-i 

o J 








a(l,0,0) 



-1,0,0) 



/ 



V 1 



o\ 




0/ 
\ 


0/ 



a(0,0, 1) = a T (0,0,-l) 



a(l,0,l) 



/ 





o\ 



















1 1 












V 1 1 





oy 








(0 





\ 






1 





1 






1 





1 











J 








( 








o\ 


-i) 


















1 
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a(0,l,l) = a T (0,-l,-l) 



a(l,0,-l) = a T (-l,0,l) 



f 








I 


1 











n 
U 


u 


u 


n 
u 










°/ 


( o 








o\ 




























1 





0/ 



a(l,-l,0) = a T (-l,l,0) 



, a(0,l,-l) = a T (0,-l,l) 



/ \ 



10 

\0 0/ 

/ \ 





\ o o i o / 

(5.3.1) 



leading to the matrix 

M(e 1 ,e 2 ,e 3 ) 

where Vj = 1 + e* 6 ^', j = 1, 2, 3. The evaluation of the determinant yields 



/ 12 -(v 2 v 3 y -( Vl v 3 y -{v x v 2 y \ 

—v 2 v 3 12 —v*v 2 -v\v 3 

—viv 3 —V1V2 12 

\ -V 1 V 2 —V1V3 —V2V3 



v^v 3 
12 J 



D{e l ,e 2 ,e 3 ) = 12 4 F{e 1 ,e 2 ,e 3 ) 

2 8 2 

F(6 U 6 2 , 6 3 ) = 1 - ~(ci + c 2 + c 3 ) - — C1C2C3 - — cic 2 c 3 (ci + c 2 + c 3 ) 

y z / oi 



+ ^(c? + c^ + C 2), 



where q = cos 2 (#j/2). This gives 

Ns T {£fcc) = 
where A = 16 3 ,n = 4iVi jV 2 iV3, so that 



Zf c 



ln(12) + 



4 



2^ 



r d9; 



2tt 



2tt 



]nF(e u e 2 ,e 3 ). 



(5.3.2) 



(5.3.3) 



nnn»(|'X'T)' (0 ' M) (5 - 3 - 4) 

fc 1= k 2 =0 k 3 =0 6 



(5.3.5) 



(5.3.6) 



The numerical evaluation of (|5.3.5|) yields the value 

z fcc = 2.354 4(4). 

6 Planar Lattices 



In this section the formulation of section 4 is applied to some other planar lattices (the square 
and 4-8-8 lattices have been considered in preceding sections). 
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6.1 Triangular Lattice 

The triangular lattice can be regarded as an Aq x N 2 square net of sites with one additional 
diagonal edge added, in the same way, to every square of the net. In this picture we have 
v = 1, a(±l, 0) = a(0, ±1) = a(l, 1) = a(-l, -1) = 1, and 

M(0i, 9 2 ) = 6 - (e ldl + e^ 1 + e i02 + e"^ 2 + e l(ei+92) + e ~ l{dl+e2) ) . (6.1.1) 

It follows that 

^ / _^i n 6-2(cos0 1 + cos0 2 + cos(0i + 2 )) 
-7T 2vr J-tt 2ix L V /J 

= ^(i_5-2 + 7 -2_ n -2 + 13 -2_ ) 

7T 

= 1.615 329 736... (6.1.2) 



z tr 



The result ( |6.1.2| ) was reported previously in where it was derived via the connection 



to the Potts partition function and a mapping to a solvable vertex model. For a finite 
triangular lattice of N\ x iV 2 cells, the number of spanning trees N$t is given by ( |4.12| ) with 
d = 2, A = v = 1, an expression we shall not reproduce here. 

6.2 Honeycomb Lattice 

It is instructive to derive z^ c using the general formulation. The honeycomb lattice is a 
square net of unit cells of u = 2 sites. Consider the honeycomb lattice in the form of a "brick 
wall" and regard the two sites connected by a vertical edges as forming a unit cell. Then 
one has 

o(0,0)=(J J), a(l,0)=a(-l,0) = a r (0,l)=a T (0,-l)= (° J), (6.2.1) 



and 



wffl ^ _ f 3 _(i + e^+e^) 

D{e u e 2 ) = 6-2fcos0 1 + cos0 2 + cos(0 1 + 2 ) N ) • (6.2.2) 



This leads to the relation Zh c = z tri /2 given in ( 1.23Q . For finite honeycomb lattices the 



number of spanning trees is the same as that of its dual, the triangular lattice. 
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6.3 Kagome and Diced Lattices 



The kagome lattice has the structure of a square net of unit cells each containing v = 3 sites 
forming a triangle. Therefore one has 



a(0,0) 



1 1 
10 1,. 

1 1 



1 





a(l,0) = a T (-l,0) =(000 

1 




a(0,l) = a T (0,-l)= , a(l, 1) = a T (-l, -1) = , (6.3.1) 



and 



M(e 1 ,e 2 ) 





4 -(l + e^ 2 ) 

(l + e iei ) -(1 + e i ^ +e '^) 



1 



with 



D(e 1} e 2 ) = i2 

This yields the result 



-(1 + e- Wl ) 
_(1 +e -K0i+6 2 y 

4 



3 — y cos 6>! + cos 6* 2 + cos(6' 1 + 9 2 

Zkag = (Ztri + ln6)/3. 



(6.3.2) 

(6.3.3) 
(6.3.4) 



For a finite kagome lattice of Ni x N 2 cells, the number of spanning trees N$t is given by 
(|JJ) with d = 2, A = 6 2 , i/ = 3. 

For the diced lattice, which is the dual of the kagome lattice, from (|1.17 ) and (|1 . 23|) , we 
obtain 

N S T(£diced) = N S T(£kag), z(C diced ) = z(£ kag ). (6.3.5) 

6.4 3-12-12 Lattice 

The 3-12-12 lattice is the lattice shown in Fig. |I|(b) which has the structure of a square net 
with unit cells each containing v — 6 sites. Label the six sites of a unit cell as shown, one 
has 



M(e l ,e 2 ) 



I 
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-1 


-1 





_ e i9 1 


o \ 




-1 


3 


-1 








_ e -i8 2 




-1 


-1 


3 


-1 
















-1 
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-1 


-1 




_ e -i9 1 








-1 


3 


-1 


\ 





_ e i9 2 





-1 


-1 


3 / 



(6.4.1) 



with 



D(9 1 , 9 2 ) = 30 



3 — ( cos 9i + cos 9 2 + cos(^i + 9 2 



(6.4.2) 
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Hence, 



3-12-12; 



1 

6 L 



z tH + ln(15) 



(6.4.3) 



For a finite 3-12-12 lattice of N\ x N2 cells, the number of spanning trees Nst is given by 
( |£g ) with A = 2 ■ 3 2 ■ 5 2 = 450 and v = 6. 

6.5 Union Jack Lattice 

The Union Jack lattice is the dual of the 4-8-8 lattice. We use ( 1.171) an d Ql-23 ) to obtain 

Nst(Cuj) = N ST (C 4 -s-8), z{Cuj) = 2 z(C A _ 8 -s), (6.5.1) 
where expressions of the right-hand sides in (|6.5.1|) are given in, respectively, ( [4.9|) and (|4.11| ). 



6.6 N x 00 Lattice Strips 

In this section we give results on the explicit evaluation of z for N x 00 strips of square and 
triangular lattices for 2 < iV < 6 by using ( |4.18| ) and ( [4.19| ). The results are 

1 



^sijr(2xoo) 

^sq(3xoo) 
z sq(4xoo) 

^q(5xoo) 
^s<jr(6xoo) 



ln(3 + 2 v / 2 ) = 0.881 373 5... 



3 V 2 
1 



1.044 532 8... 



ln(3 + 2V2 ) + 2 ln(2 + V3 ) 
9 - + (70 - 18V$) 1/2 



In 



4 
2 
5 

1.123 728 0... 
1 



1.099 165 7... 

'9 + v / 5 + (70 + 18V/5) 172 



In 



4 



2„ 1( i±^Lu2i„r 5 



6 



1.136 865 4... 



+ ln(3 + 2V2 ) 



(6.6.1) 



and 



^tri(2xoo) 
^tri(3xoo) 

^tri(4xoo) 
^tri(5xoo) 



21n2 = 1.386 294 3.. 
1 



-ln2 + 21n(7 + 3v / 5 ) = 1.514 280 5... 
21n2 + ln(3 + ^7 ) = 1.558 598 8... 
-7hi2 + 21n(l3- v / 5 + (150 - 34^ ) 1/2 
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+2 lnfl3 + Vh + (150 + 34^ ) 1/2 



1.579 041 2. 



ln(5 + Vl~3 ) + ln(7 + 3\/5 



1.590 133 9.. 



(6.6.2) 



One observes that the values of z are monotonically increasing in N, in accordance with The- 
orem 2.1. One also observes that zc(n 2 xoo) converge reasonably quickly toward the respective 
values z sq in (|2.17j) and z tri in ( |2.19| ). For example, z sq (^ xoo ) and z sq (Q XOO ) are, respectively, 
within 10 % and 2.5 % of the value ( p.!7| ) for the infinite square lattice. Similarly, Zt r i(3x<x>) 
and £4^(6x00) are, respectively, within 6 % and 1.5 % of the value ( |2.19| ). 



6.7 Connection with the Critical Ising Model and Dimers 

In this section we establish a result relating z^ for planar lattices £ to the free energy of the 
Ising model on C at the critical point. We also remark on a connection of spanning trees 
with dimers for planar lattices. We first state our result as a theorem. 



Theorem 6.1. 

For planar lattices C we have the identity 

zc = ac + Vl (6.7.1) 

where ac is a lattice-dependent constant, and is the (reduced) free energy of the Ising 
model on £ at the critical point, 

Proof. We use the fact that the Ising model is the infinite bare quartic coupling (A — > 00) 
limit of the 4 lattice field theory fl2~5 |. The partition function for the <fr quantum field 



theory is 

Z= I |n^| e " 5 ( 6 - 7 - 2 ) 



where the action S is an integral of the Lagrangian density with the quadratic part 



2 



^-4/ A £(f) 2+m2 4 (6 - 7 - 3) 

j=i 3 

After integrating by parts, the kinetic terms becomes (1/2) J d 2 x(f)[—d 2 + m 2 ]0, where d 2 
is the Laplacian. Further discretizing to a lattice £, the integrand in (|6.7.3| ) becomes the 
summand 
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where Qij are elements of the Laplacian matrix Q given by (JO 

At the critical point, because the phase transition is of second order with a divergent 
correlation length, the mass m in ( |6.7.4|) which is the inverse correlation length vanishes. 
One is left simply with the term involving Q. Letting A — > oo, the functional integrals 
are now reduced to discrete sums over the Ising variables Oi = ±1, and, from the Onsager 
solution and the correspondence with (|4.15|) , one finds 

d6 l r d6 2 



fl 



2 2v 



2tt 



2tt 



ln£>( 



(6.7.5) 



where ac is a lattice-dependent constant. This establishes (|6.7.1| ) after using ( |4.16| ). □ 



We find the explicit results 



Zsq 

Zhc 
Zkag 



-ln2 + 2/ s c g 
ln(V3/2) + 1ft n 
- ln(2^3 ) + 2f c hc 

Iln3-~ln2-~ln(2 + ^3)+2/, 



kag 



(6.7.6) 



The values of in (|6.7.6| ) are well-known |26|] 



Finally, we remark that Temperley |27]] has established a bijection between spanning trees 
on an N x N square lattice with free boundaries and dimer configurations on a (2iV — 1) x 
(2iV — 1) square lattice with one boundary site removed. We shall not repeat the proof here, 
except to note that this equivalence can be extended more generally to arbitrary planar 
graphs p8 , |29|. Together with Theorem 6.1, this bijection implies a connection between 
dimers and the critical Ising model, a connection which has been observed by Fisher in a 
weaker form . 



7 Discussion 



It is of interest to investigate how close the actual value of zc is to the Mckay-Fan-Yau upper 
bound ( |2.7|) . To do this, we define for K-regular lattices the ratio 



lnC K 



(7.1) 



where C K is given by ( [2.6|) . For lattices which are not K-regular, we compare zc instead with 
the general upper bound ( |2.4|) and consider the ratio 

zc 



rc 



In k 



eff 



(7.2) 
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Table 1: Values of zc and rc for different lattices C. 



c 


d 




K eff 




zc 


rc 


3-12-12 


2 


3 




1/2 


0.720 563 3... 


0.861 


4-8-8 (Bathroom-tile) 


2 


3 




1/2 


0.786 684(1) 


0.940 


Honeycomb 


2 


3 




1/2 


0.807 664 8... 


0.965 


Kagome 


2 


4 




1 


1.135 696 4... 


0.933 


Diced 


2 




4 


1 


1.135 696 4... 


0.819 


£2 (Square) 


2 


4 




1 


1.166 243 6... 


0.959 


Union Jack 


2 




6 


2 


1.573 368(2) 


0.878 


Triangular 


2 


6 




2 


1.615 329 7... 


0.955 


£3 (Simple cubic) 


3 


6 




2 


1.674 148 1(1) 


0.990 


Cbcc (Body-centered cubic) 


3 


8 




3 


1.990 2(1) 


0.991 


u 


4 


8 




3 


2.000 0(5) 


0.996 


£5 


5 


10 




4 


2.243 


0.998 




6 


12 




5 


2.437 


0.999 


Cf cc (Face-centered cubic) 


3 


12 




5 


2.354 4(4) 


0.965 


£focc(4) 


4 


16 




7 


2.732 3(1) 


0.998 



Results are summarized in Table |I[ 

As is evident from Table 1, for regular lattices that we have studied, zc is a monotonically 
increasing function of the coordination number n. We also observe that, for a fixed k, the 
value of zc increases with the spatial dimension of the lattice zc- Examples are (i) the 
triangular and simple cubic lattices (both with k — 6), (ii) the bcc and d = 4 hypercubic 
lattices (k = 8), although the difference between Zb cc and zc 4 is very small, and (iii) the 
fee and d = 6 hypercubic lattices (n = 12). Furthermore, our results indicates that in two 
dimensions the square lattice is a little more densely connected than the kagome lattice, 
both of which have k = 4. The ratio rc is observed in several cases to increase with k, but 
not in all cases; a counterexample is r = 0.959 for the square lattice (k = 4) and r = 0.955 
for the triangular lattice (re = 6). 
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